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Weak values allow the measurement of observables associated with noncommuting operators. Up 
to now, position-momentum weak values have been mainly developed for (relativistic) photons. 

In this Letter, a proposal for the measurement of such weak values in typical electronic devices 
is presented. Inspired by the Ramo-Shockley-Pellegrini theorem that provides a relation between 
current and electron velocity, it is shown that the displacement current measured in multiterminal 
configurations can provide either a weak measurement of the momentum or strong measurement of 
position. This proposal opens new opportunities for fundamental and applied physics with state-of- 
the-art electronic technology. As an example, a setup for the measurement of the Bohmian velocity 
of (nonrelativistic) electrons is presented and tested with numerical experiments. 


Introduction. —Nowadays, there is a rapidly growing 
interest in weak measurements and weak values [1-4], 
both from fundamental and applied points of view. Since 
weak values (a weak measurement postselected by a 
strong measurement) provide information on incompati¬ 
ble observables associated with noncommuting operators, 
relevant topics of quantum mechanics, such as the tun¬ 
neling times [5], Hardy’s paradox [6, 7], Leggett-Garg in¬ 
equalities [8, 9], and quantum amplification [10-12], are 
being revisited. Especially attractive is the simultane¬ 
ous measurement of position and momentum: a set of 
weak measurements of position postselected by a strong 
measurement of momentum is proportional to the wave 
function of the system [13, 14], while a weak measurement 
of momentum postselected by a strong measurement of 
position gives the local velocity of a quantum particle 
[15, 16]. 

Most experimental techniques for weak values are de¬ 
veloped for photons, whose technology is not easily trans¬ 
ferable to industry based on electronics. The few propos¬ 
als dealing with weak measurements in solid-state sys¬ 
tems [17-22] use particle current measurement (i.e., elec¬ 
tron charge detection). Instead, we propose measuring 
displacement current (i.e. time dependent variations of 
the electric field) to get information on the position and 
momentum of a quantum state. Similar to Landauer’s 
proposal [23] which demonstrates that the measured dc 
current provides information of the quantum transmis¬ 
sion coefficient, here, we show that the weak measure¬ 
ment of the ac current flowing through a (properly pre¬ 
pared multiterminal) electron device provides informa¬ 
tion on the whole quantum state. This new proposal 
opens original routes to study, both, fundamental physics 
and quantum engineering. 

As an example of the potentialities of our proposal, 
inspired by the old classical works of Shockley and 
Ramo [24, 25], we discuss the measurement of the local 
(Bohmian) velocity (i.e. the current density divided by 
the modulus of the wave function) for an electron. Such 
velocity is obtained from a weak value constructed from 
two measurements of the displacement current on two 
different metallic surfaces belonging to a multiterminal 


device. The electric field generated by a moving electron 
(which contains information on the electron dynamics) is 
detected in a large metallic surface even when the elec¬ 
tron is far from that surface. Such metallic surface (i.e. 
the electrons inside) does only weakly perturb the quan¬ 
tum state of the electron. A strong measurement of posi¬ 
tion can be envisioned by using a small surface that only 
detects the electric field when the electron is very close. 
Next, before describing the simpler strong measurement 
for postselection, we explain the weak measurement. 

Weak measurement of the total current .—The mea¬ 
surement of the electrical current can be understood as 
a two step process. The first step is an electromagnetic 
propagation of the total current along the cable (that 
connects the quantum system and the ammeter in Fig. 
1). The total current through Si is equal to the current 
through Sa far from the active region. This equivalence 
(due to the divergencelessness of the total current) is ex¬ 
act for the sum of the particle plus the displacement cur¬ 
rents, but not for particle current alone. The second step 
is done by the ammeter that transforms the total current 
Sa into a pointer value. 

In this Letter we analyze this two-step many-body 
quantum measurement with the corresponding quantum 
errors and backaction [26]. We define the density matri¬ 
ces of the system and the probe, at the initial time to, as 
Psys{to) and Pproifo) 1 respectively. The time evolution of 
the density matrix of the entangled system between the 
initial time to and the final time t m is: 

Ptotifm ) = U(t 0 ,t m ) (p{t-o)sys & p(fo)pro) U (to, t m )^l) 

where the unitary operator U(to, t m ), which contains the 
free part of the system and of the probe as well as the 
electromagnetic interaction between them, is the respon¬ 
sible of the first step, i.e., translating the total current 
from Si towards Sa (see Fig. 1). 

The second step of the measurement is done at time 
t m by a projective measurement of the probe performed 
by the ammeter. Such measurement provides the value I 
on the pointer [27]. Thus, the probability of the specific 
result I associated with the eigenstate | J) in the ammeter 
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at time t m is: 

= Tr pro (\I){I\p' pro {tm)) , (2) 

where p' pro (t m ) = Tr ays (p tot (t m )) is the reduced density 
matrix of the probe, with Tr sps the partial trace opera¬ 
tion over the system coordinates. Hereafter, we will con¬ 
sider only one particle in the system, while considering 
an arbitrary large number Np of particles in the probe. 
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FIG. 1. Schematic representation of the whole setup divided 
into three regions. The system is enclosed in the volume Q 
(red dashed lines) with an electron inside. The left and right 
probes (cables) with Np electrons ensuring that the total cur¬ 
rent on the left surface Si is equal to that on the left Sa 
(green dashed-dotted lines). Finally, the left and right meters 
(ammeters) indicate the measured value of the total current 
(the right probe and meter will be substituted in Fig. 3(b) by 
a multiterminal structure). 

The standard way of describing the probability 
V(I,t m ) is not by referring to the whole probe’s coor¬ 
dinates in Eq. (2), but only to the quantum system’s co¬ 
ordinates: 

nl tm) = Tr sps (l{i)p' sys {t m j) , (3) 

where p' sys (tm) = Tr pro ( / 5 tot (t m )) is the system reduced 
density matrix, with Tr pro the partial trace operation 
over the probe’s coordinate. In Eq. (3), we define 1(1) as 
a general positive operator-valued measure (POVM): 

I(I) = j \I)g(I,I)(I\dI, (4) 

with g(I,I) a real positive number. If we intro¬ 
duce the quantum system density matrix p' sys (t m ) = 
J2jPj(.t m )\i>j)(ipj\ with pj(t m ) a normalized probability 
and \i/jj) = f aj(t m )\l)dl into Eq. (3), together with Eq. 
(4), the term g(I , I) can be computed from the relation: 

P(j,tm)=5>(*) / 9(i,I)\a J i(tm)\ 2 dI. (5) 

J 

Next, we discuss in what conditions the measurement of 
the total current provides information on the momentum 
of the quantum system. 

Displacement current and momentum measurement .— 
The (quantum ensemble) value of the total current (I) 


can be computed [28, 29] straightforwardly as the sum of 
particle plus displacement current on the surface Si: 

(I(t))= [ (Jc( r , t)) ■ ds + [ e d(E ^ )} - da, (6) 

J Si j Si 

where e is the dielectric constant, (J c (r,f)) the quantum 
ensemble value of the quantum particle current density 
and (E(r, £)) the quantum ensemble value of the elec¬ 
tric field on the surface. Identically, from the Ramo- 
Schockley-Pellegrini theorem with a quasistatic approxi¬ 
mation that neglects the vector potential (magnetic) con¬ 
tribution [28-31], such a mean value in Eq. (6) can be 
written as 

(I(t))=-f F(r)-(J c (r, t))-dv+f e-F(r)- *)) , rfs ,(7) 
aJ n J s 

where S is the close surface of the volume f l (see Fig. 
1) that contains Si and (V(r, £)) is the ensemble value 
of the electrostatic potential. In the particular situation 
where electron transport takes place between two metal¬ 
lic surfaces included into O, no variation of the potential 
appears on these surfaces, i.e., d(V(r,t))/dt = 0. If the 
square of the distance between the metallic surfaces is 
smaller than the area of the surfaces, L 2 -C Si, we get 
F(r) x = 1 /L x [28, 32]. Then, Eq. (7) is rewritten as 

m) =-rj-mh, ( 8 ) 

TflL/x 

where (p(£))n is the mean value of the momentum in 
the volume D of the device. We have assumed that the 
support of the density matrix of the quantum system is 
inside the volume ft. See Ref. [26] for an alternative 
demonstration of Eq. (8). The experimental evaluation 
of this mean momentum requires the knowledge of the 
probabilities in Eq. (5). 

Numerical evaluation of V (/, t). —In order to compute 
the probability V(I,t), one needs to simulate the time 
evolution in Eq. (1) and use Eq. (2). Obviously, the exact 
solution of this problem is not accessible, but fortunately 
we can provide some reasonable simplifications to handle 
this problem. The many-particle Coulomb interaction 
among electrons is enough to consider the transmission 
of the total current from Si towards Sa [29, 30, 33]. 

To numerically treat the problem we use the 
method reported in Ref. [34] based on the def¬ 
inition of conditional wave functions if?(xi,t) = 
4t(X{(£), ...,X J i _ 1 (t),Xi,X 3 i+1 (t), ...,£) for each t-particle 
with Ik the many particle wave function [34-37] . For sim¬ 
plicity, the quantum system is treated as a ID system. 
The capital letter X^(i) denotes the actual (Bohmian) 
particle 3D position. The subindex i = 1,iVp + 1 refers 
to the particle in the quantum system plus the Np par¬ 
ticles in the probe. The superindex j = 1 denotes 

one (of the infinite) particular selections of the initial po¬ 
sitions. The time evolution of ipi(xi ,£) is obtained by 
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solving numerically the following single-particle (condi¬ 
tional) Schrodinger-type equation: 

= [.Hq + V] V’j(aq, t) (9) 


where V = V(xi, X^f), X. 3 Np+1 (t)) is the conditional 
Coulomb potential felt by the system and Hq is its free 
Hamiltonian. The trajectory X[{t) is obtained from 
the so-called guidance equation vf(t) = dX{{t)/dt = 
(fi/TO)Im((<9^/clzq)(l/'f/>()). This numerical method 
[34] for our problem allows the following: (i) A manage¬ 
able treatment of the many-particle interaction of Eq. (1) 
for a small interaction [38] and (ii) a simple treatment of 
Eq. (2) by using, for each experiment, one channelized 
conditional wave function \&(aq, X^t),..., X^- +1 (t), t) 
[37]. The quantum probabilities computed from Bohmian 
trajectories are, by construction, identical to the ones ob¬ 
tained from standard quantum tools. Once an (infinite) 
j = 1 ensemble of initial positions are considered, 

we compute the probability V{I , t) in the range [7, 7+77] 
as V(I, t)dl = i - 7)0(7 + 77- P{t)), with 

O(x) the Heaviside function. The measured total current 
P(t ) is computed as: 


P(t) 



1+Np 

■ds= £ eV$(X J fc (t))-v J fc (t),(10) 

k =1 


where E J is the electric field generated on the surface Si 
by the 1 + Np particles X. J k . We include the previous 
conditional wave function algorithm in the right-hand 
side of Eq. (10) where $ is the electric flux due to the 7th 
electron and v 2 ( t ) is the local (Bohmian) velocity. We 
shall use Eq. (10) as a numerical tool for calculating the 
subsequent results (see details in Ref. [26]). 

In Fig. 2 we show the numerical results of V(I,t m ) 
when M = 55000 for a large surface Sa, he., for Sa 
L 2 . For simplicity, the electron 1 in the quantum system 
is moved only in the ID transport direction x through 
Eq. (9) interacting, through Coulomb potentials, with 
the other electrons in the left and right probes. The 
Np electrons are simulated semiclassically in a 3D space 
taking into account the many-body Coulomb interaction 
among them (and with the particle in the system) plus 
the interaction with a bath of phonons at room temper¬ 
ature. Although we obtain a large dispersion of values of 
7, its mean value exactly coincides with the mean value 
(7) obtained without including the Np electrons in the 
metals. In addition, we observe that the system is only 
slightly modified by the interaction with the electrons in 
the metal when the distance is reduced (see inset in Fig. 
2 )- 

We use Eq. (5) to find g(I,I). The support of the 
function V (7, t) is much larger than the current (mo¬ 
mentum) distribution of the wave function, |a/(f)| 2 for 
any J. This information can be obtained comparing the 
distribution V(I,t) in Fig. 2 and the |o/(t)| 2 , which are 



I(M) 


FIG. 2. Probability distribution (red solid line) of the mea¬ 
sured total current from Eq. (10) from M = 55000 numerical 
experiments. Mean value of the distribution (green dashed 
line) obtained from the same simulations without the mea¬ 
surement backaction. Inset: Error on the wave function com¬ 
puted as Erro r wave = J \ip(xi, t m ) — ^^{xi^m^dxi at 
the final time t m for different values of the distance d in Fig. 
1. The wave function i/j(xi,t m ) is computed from Eq. (9) with 
the many-particle potential, while ip rnearl (xi,tm) from Eq. (9) 
with only an external mean-field potential. 


known directly from the wave function \ipj). Therefore, 
we can approximate |a/(t)| 2 ~ S((I) — I) in Eq. (5) to get 
g(I, (7)) « V(I,t m ). We have used = 1 - From 

these last results, along with Fig. 2 and Eq. (8), we obtain 
a specific expression for the ID version of the operator 
in Eq. (4): 

Lipw) « C w Jdpe 2 < \p)(p\, (11) 

where a w is the width of the Gaussian distribution in Fig. 
2 (related to the specific plasma oscillations of the metal 
probes) and C w is a suitable normalization constant. The 
subindex w reflects the particular geometry of the device 
that leads to a weak measurement. 

Measurement of Bohmian velocities. —Hereafter, to 
show the possibilities of these weak measurements of the 
displacement current, we provide a proposal for measur¬ 
ing the Bohmian velocity of a quantum particle. We con¬ 
sider a multiterminal device where there are two metal¬ 
lic surfaces working as sensing electrodes. See Fig. 3(b). 
The length of the device is L x = 280 nm and the sur¬ 
face S w « 10“ 11 m 2 , satisfying the relation S w L 2 . 
The total current can be obtained from the operator in 
Eq. (11). The other electrode is divided into n surfaces, 
each one electrically insulated from the others, connected 
to its s—ammeter and satisfying the opposite relation 
S Si <C A 2 . In these small surfaces the total current in Eq. 
(7) can be computed from F(r) x = — ae~ a ^ r ~ rs ^ with 
a = y/2 /S Si and r s = {a: s , y s , z s } the central position 
of S 8i [32]. Here, Eq. (7) provides a measurement of the 
total current that becomes different from zero when the 
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particle is going to cross one of the small surfaces of Fig. 
3(b). In this small surface, the ID version of the total 
current operator in Eq. (4) is: 

r (*-3- ? ) 2 

I s {x s ) « C s dxe 2<t - \x)(x\. (12) 


where is related to S Si in Fig. 3(b). See Ref. [26] for 
an alternative derivation of Eq. (12). 

The exact procedure for measuring the Bohmian veloc¬ 
ity is the following: The total density matrix evolves from 
to till t m following Eq. (1). At time t m the lu-ammeter 
measures weakly the momentum of the particle through 
Eq. (3). Then, the quantum system evolves until a peak 
of current is measured (or not) in the Si-ammeter. The 
values measured with the w-ammeter will be postselected 
by the measurement of the positions performed by the Si- 
ammeter. Thus, we compute 


E[p u 


f dp w p w V(p w n x s ) 
V{x s ) 


(13) 


which under the condition a w ( h/a s ) becomes 


E \p w \x s \ _ J(x s ,t m ) = v ( x ^ tm) (14) 


\lt>(x s ,tm )\ 2 


where m is the particle mass and v(x s , t m ) is exactly the 
Bohmian velocity. See Ref. [26] for an explicit develop¬ 
ment from Eq. (13) till Eq. (14). 

Hereafter, to provide a realistic estimation on how 
many experiments M are needed to capture typical quan¬ 
tum interference phenomena in electronic devices and to 
understand the backaction in the system, we provide a di¬ 
rect numerical simulation of the whole weak value proce¬ 
dure by simulating Eq. (1) from to till t m and a posterior 
measurement through Eq. (3) using the conditional wave 
function technique mentioned in Eq. (9), without using 
the operators defined in Eq. (11) and Eq. (12). In Fig. 
3(a) it is reported the Bohmian velocity v(x s ,t ) obtained 
from an ensemble of M = 55000 (identically prepared) 
two-time measurement experiments calculated from Eq. 
(14). The different errors are due to the different num¬ 
ber of particles that effectively arrive at each position. 
In order to see interference effects we consider as initial 
wave function ^>(a;i,0) a superposition of two Gaussian 
wave packets, in the motion direction x, whose central 
positions are separated of 50 nm at the initial time (see 
Fig. 3c). Each Gaussian wave packet has the same dis¬ 
persion of 3 nm and energy of 0.0905 eV. The velocity 
field exhibits the typical interference pattern. 

In Fig. 3(c) we report the comparison of the trajecto¬ 
ries obtained from the described procedure and the wave 
function of the single particle problem. As expected, the 
trajectories are more dense near the maximum and less 
dense near the minimum of the interference pattern. 

Conclusions. —In conclusion, we have presented an ex¬ 
perimental proposal for measuring weak values of posi¬ 
tion and momentum in solid state devices. In particular, 
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FIG. 3. (a) Bohmian velocity (red points with errors) from 
Eq. (14) as a function of position x s at time t m = 0.3 ps, ob¬ 
tained from an ensemble of numerical experiments with the 
backaction of the ammeter. Its errors are calculated from 
the standard deviation of the set of velocities at time t m and 
position x s divided by the square root of the number of trajec¬ 
tories passing through position x s at time f m . The Bohmian 
velocity (green solid line) in position x a obtained from the 
same set of simulations without the backaction. (b) Schematic 
representation of the multiterminal device described in the 
text where S w and different S Si are indicated, (c) Wave func¬ 
tion (upper plot) obtained from simulations without the back- 
action of the measurement and Bohmian trajectories (lower 
plot) reconstructed from the M experiments with backaction, 
at different times. The green dashed line represents the time 
chosen for the plot in a). 


we show how the Bohmian velocity of a (massive) elec¬ 
tron can be obtained from weak values of the displace¬ 
ment currents in a multiterminal electron device. We 
emphasize that the Wiseman’s protocol [15] for (weak) 
measuring the Bohmian velocity was developed for mas¬ 
sive nonrelativistic particles (not for massless relativistic 
ones). So, our Bohmian velocity measurement for elec¬ 
trons in solid-state structures exactly fulfills the nonrela¬ 
tivistic scenario contemplated in Ref. [15], while the local 
velocity of an ensemble of relativistic photons measured 
in the pioneering work of Kocsis et al. [39] does not. The 
feasibility of our proposal has been tested numerically in 
the iconic double-slit experiment at room temperature. 
In real experiments, for samples at low temperatures with 
few microns coherence length [40] , the protocol presented 
here can be implemented with frequencies below 50 GHz. 
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Our protocol implies the use of the (now available) single 
electron sources [41, 42], while standard thermal injec¬ 
tion in solid-state devices would provide the (ensemble) 
velocity of the mixed state [43]. This work opens the 
path for answering intricate fundamental questions and 
developing new quantum engineering applications using 
the successful electronic/semiconductor industry. 

The authors would like to thank Tomas Gonzalez, 
Javier Mateos, Philippe Dollfus and Massimo Macucci for 
fruitful discussions and suggestions. This work has been 
partially supported by the Fondo Europeo de Desarrollo 
Regional (FEDER) and Ministerio de Economa y Com- 
petitividad through the Spanish Projects No. TEC2012- 
31330 and No. TEC2015-67462-C2-1-R, the General- 
itat de Catalunya (2014 SGR-384), and by the Eu¬ 
ropean Union Seventh Framework Program under the 
Grant Agreement No. 604391 of the Flagship initiative 
“Graphene-Based Revolutions in ICT and Beyond”. N. 
Z. is supported in part by INFN. 


* xavier.oriols@uab.cat 

[1] Y. Aharonov, D. Z. Albert, and L. Vaidman, Phys. Rev. 
Lett. 60, 1351 (1988). 

[2] A. G. Kofman, S. Ashkab, and F. Nori, Phys. Rep. 520, 
43 (2012). 

[3] B. Braverman and C. Simon, Phys. Rev. Lett. 110, 
060406 (2013). 

[4] J. Dressel, M. Malik, F. M. Miatto, A. N. Jordan, and 
R. W. Boyd, Rev. Mod. Phys. 86, 307 (2014). 

[5] N. Vona, G. Hinrichs, and D. Diirr, Phys. Rev. Lett. Ill, 
220404 (2013). 

[6] J. S. Lundeen and A. M. Steinberg, Phys. Rev. Lett. 102, 
020404 (2009). 

[7] K. Yokota, T. Yamamoto, M. Koashi, and N. Imoto, New 
J. Phys. 11, 033011 (2009). 

[8] A. Palacios-Laloy, F. Mallet, F. Nguyen, P. Bertet, D. 
Vion, D. Esteve, and A. N. Korotkov, Nature Phys. 6, 
442 (2010). 

[9] J. Dressel, C. J. Broadbent, J. C. Howell, and A. N. 
Jordan, Phys. Rev. Lett. 106, 040402 (2011) 

[10] O. Hosten and P. Kwait, Science 319, 787 (2008). 

[11] P. B. Dixon, D. J. Starling, A. N. Jordan, and J. C. 
Howell, Phys. Rev. Lett. 102, 173601 (2009). 

[12] D. J. Starling, P. B. Dixon, A. N. Jordan, and J. C. 
Howell, Phys. Rev. A 82, 063822 (2010). 

[13] J. S. Lundeen, B. Sutherland, A. Patel, C. Stewart, and 
C. Bamber, Nature 474, 188 (2011). 

[14] J. S. Lundeen and C. Bamber, Phys. Rev. Lett. 108, 
070402 (2012). 

[15] H. M. Wiseman, New J. Phys. 9, 165 (2007). 

[16] F. L. Traversa, G. Albareda, M. Di Ventra, and X. Oriols, 
Phys. Rev. A 87, 052124 (2013). 

[17] A. Romito, Y. Gefen, and Y. M. Blanter, Phys. Rev. 
Lett. 100, 056801 (2008). 

[18] N. S. Williams and A. N. Jordan, Phys. Rev. Lett. 100, 
026804 (2008). 

[19] V. Shpitalnik, Y. Gefen, and A. Romito, Phys. Rev. Lett. 
101, 226802 (2008). 

[20] O. Zilberberg, A. Romito, and Y. Gefen, Phys. Rev. Lett. 


106, 080405 (2011). 

[21] A. Bednorz and W. Belzig, Phys. Rev. Lett. 105, 106803 

( 2010 ). 

[22] C. Laflamme and A. A. Clerk, Phys. Rev. Lett. 109, 
123602 (2012). 

[23] R. Landauer, IBM J. Res. Dev. 1, 223 (1957) 

[24] W. Shockley, J. App. Phys. 9, 635 (1938). 

[25] S. Ramo, Proceedings of the I. R. E. 27, 584 (1939). 

[26] See the Supplemental Material [url] for a discussion of 
technical details relevant to the main text. 

[27] We emphasize that the ammeter, as any other electronic 
devices, works as a low-pass filter. Therefore, the current 
measured by the ammeter, after the electronic filtering, 
will be Ifreq(t ) = 1/T f* +T I(t')dt', i.e. the current mea¬ 
sured at frequency / = 1/T. Our discussion is valid for 
frequencies / > 1/r, where r = L x /v is the dwell time of 
the electron in the device, which depends on the length 
of the device L x and on the velocity v of the electron. 
In the text we have defined I = If req (t) to simplify the 
notation. 

[28] B. Pellegrini, Phys. Rev. B 34, 5921 (1986). 

[29] G. Albareda, F. L. Traversa, A. Benali and X. Oriols, 
Fluct. Noise Lett. 11(3), 1242008 (2012). 

[30] The spatial dependence of the vector potential can be 
neglected in front of the spatial dependence of the scalar 
potential for active regions of nanometric dimensions un¬ 
til frequencies of few hundreds of THz. 

[31] M. Lanzagorta, Quantum Radar (Morgan & Claypool, 
San Rafael, 2011). 

[32] A. Benali, F. L. Traversa, G. Albareda, A. Alarcon, 
M. Aghoutane and X. Oriols, Fluct. Noise Lett. 11(3), 
1241002 (2012). 

[33] In particular, we will concentrate only in the interac¬ 
tion between the electron of the quantum system and 
Np electrons in the probe which are close enough to the 
system (the coulomb force is inversely proportional to the 
squared distance). 

[34] X. Oriols, Phys. Rev. Lett. 98, 066803 (2007). 

[35] D. Diirr, S. Goldstein, and N. Zanghi, J. Stat. Phys. 116, 
959 (2004). 

[36] D. Diirr and S. Teufel, Bohmian Mechanics (Springer, 
2007). 

[37] D. Diirr, S. Goldstein, and N. Zanghi, J. Stat. Phys. 67, 
843 (1992). 

[38] G. Albareda, J. Sune, and X. Oriols, Phys. Rev. B 79, 
075315 (2009). 

[39] S. Kocsis, B. Braverman, S. Ravets, M. J. Stevens, R. P. 
Mirin, L. K. Shalm, and A. M. Steinberg, Science 332, 
1170 (2011). 

[40] F. Miao, S. Wijeratne, Y. Zhang, U. C. Coskun, W. Bao, 
and C. N. Lau, Science 317 1530 (2007). 

[41] E. Bocquillon, V. Freulon, J.-M Berroir, P. Degiovanni, 
B. Pla§ais, A. Cavanna, Y. Jin, and G. Feve, Science 339 
1054 (2013). 

[42] J. Dubois, T. Jullien, F. Portier, P. Roche, A. Cavanna, 
Y. Jin, W. Wegscheider, P. Roulleau, and D. C. Glattli, 
Nature 502, 659 (2013). 

[43] A. Luis and A. S. Sanz, Ann. Phys. 357, 95 (2015). 



6 


A - DERIVATION OF THE BOHMIAN 
VELOCITY FROM KRAUS OPERATORS 


it is possible to work out the integral between squared 
parenthesis in Eq. (18), 


Alternatively to the numerical experiment developed in 
the letter (see Fig. 3 in the letter), hereafter we provide a 
detailed analytical derivation of Eq. (14) using the Kraus 
operators I w and I s defined in the letter, Eq. (11) and 
Eq. (12), respectively. We rewrite the final expression 
here as: 


E\p w \x s \ 1 f dp w p w V(p w r\x s ) 
v{x s ,t m ) = -=-—r-• (15) 


m 


m 


V(x s ) 


First, we calculate the probability V(x s ) in denomina¬ 
tor of the last expression as: 


V{x s ) = J dp w V(p w C\x s ) 

= f dp w dp w WiuLJtfsUtJu 


dp'dp" 


I dp w 6 


. (p'-Pm) 2 _ (p"-P vj) 2 


■ wwiu'PJsUtjjf'W' 


(16) 


where in the last line it has been used the definition in 
Eq. (11) and Eq. (12). We can focus only in the integral 
between squared parenthesis in the last line of Eq. (16) 


f _ ip'-Puj ) 2 _ (p"-P-w ) 2 

Cl / dp w e 2 °l e 2 "» = 


— C'w I dp w e 


-+PM -zr + 


= e 


(p'-p") 2 


(17) 


where it has been used C w = (yCra w ) 1 / 2 . Thus Eq. 
(16) becomes: 


Ci / dxe 




h 1 h ) 1 ‘ 2mti~ 2 mh 


= e -5fr(p'-A') 2 e -|x sP '+i^r e i^p"-i^ = 

= e _ ^ (p_p ) (p'\Uf m \x s }(x s \Ut m \p")- 

Therefore Eq. (18) becomes: 


(19) 


r r (. p ’- p ") 2 , , „, 2 

V(x s ) = dp dp"e e p ) • 

• (i>\p , )(p , \u} m \x s }(x s \u tm \p")(p"\^), (20) 

if the condition 


h 

&W •) 

Or 


( 21 ) 


is satisfied then Eq. (20) simply becomes: 


V(x s ) = (V>| = mx s ,t m )\ 2 . (22) 

The second step is the calculation of the numerator of 
Eq. (15): 


j dp w p w V{p w n x s ) — 

= JJ dp'dp" Cl J dp w p w e 
■ W){p'\U}Jli s U tm \p"){p"\^), (23) 


{p' —Pvj) 2 _ (p" ~Pw ) 2 1 

2cr^ g 2cr^ 


again it is possible to calculate the integral between 
squared parenthesis in the last expression, 


V(x a ) = 


. ( P' —p") 


dp'dp" e (iji\p')(p'\Ui PJ s U tm \p")(p"\ip) 


dp'dp" (ip\p')e ^ (p"\ 


■ C 2 S J dxe {p'\Uljx){x\U trn \p") 

Recalling that: 

• (i) U tm = J dp\p)(p\e~ l! ^it, 


• (ii) (x\p) = J— ~ e lP Z 


/hvh 


(18) 


(iii) (x\UtJp") = {x\Jdp\p){p\e 1 wf | p ") = 


1 


p +i p_^_ 


/2irh 


e 6 2 mh 


r _ (p'-P vj ) 2 _ (p"-P vj ) 2 

Cl / dp w p w e 2<r » e 2<t ™ 


_ (p'-p'') 2 r 

= e ^ Cl / dp w p w e 


p ' +p »^ e _iPCjpi 


. (p'+p") 1 


Then Eq. (23) becomes: 


J^ dp w p w v(p w n x s ) — 


D r jj" \ _ (p'-p") 2 

dp'dp" [ 1 e 


W\ P ')(p'\ulitiA m \p")(p"W, 


(24) 


(25) 
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making the same steps done from Eq. (18) till Eq. (20), 
Eq. (25) becomes: 


J d'Pw'Pw^PiP'u: 11 Xg) — 
p' + p" 


= / / dp'dp" 


(p'-p ") 2 p 2 2 


• W’\p , )(p , \u} m \x s }{x s \Ut m \p")(p"\ijj). (26) 

If Eq. (21) is satisfied then Eq. (26) becomes: 


(9) in the text through the use of the Conditional wave 
function algorithm described in the letter. 

According to our discussion in the letter, the cap¬ 
ital letters {Xf(t),Y?(t),Z{(t)} denotes the actual 
(Bohmian) positions of the particles, where i identifies 
the * — th particle and j define one (from the infinite) 
particular selections of the initial positions. The flux of 
the electric field through a general ideal surface S w (see 
Fig. 1), defined as a plane of area L y ■ L z perpendicular 
to the x direction and placed in x = x w , i.e. defined by 
the points {2^,0 < y' < L y , 0 < z' < L z }, generated by 
a particle in position {X , Y, Z} can be calculated as: 


dp w p w V(p w n x s ) = 

'mu]Jti s u tm \i>) + mlJlisUtM) = 

= Re (&\U}Jtl s U tm m) (27) 

where in the last expression it has been used the prop¬ 
erty p = f p\p)(p\dp. If one writes the momentum op¬ 
erator in position representation easily realizes that the 
expression between parenthesis in the last line of Eq. (27) 
is 


Q(X,Y,Z)= I V(X,Y,Z,x w ,y',z')-ds , (30) 

Js w 

where we have eliminated the superindex j, subindex i 
and time t to simplify the notation. The electric field E 
is just computed from the Coulomb force of the electron 
in the mentioned surface. In the simple case in which the 
particle is located in {X,L y /2,L z /2}, and it moves only 
in the x direction and L y = L z = L (S w = L 2 ), Eq. (30) 
becomes: 


(ilj\U^ m ili s U trri p\ip) = -ih'i/;*(x s ,t m )-Q^4’(x s ,t m ) 


(28) 


and thus its real part is 


Re ((MU} lUA m P 


't m ± 8 ± S 

H f dlp(x s ,t m ) 


2i \ dx s 
= mJ(x s , t m ). 


1p*(x s ,tm) - (x s ,t m ) ^ x ^t m ) 


dx s 


(29) 


Therefore, using Eq. (22) and Eq. (29), the right hand 
side of Eq. (15) becomes: 


<h(X) = —tan 1 [ - w - . 

\A{x w -X)^{x w -Xy + ?f] 

(31) 

Let us evaluate Eq. (31) in the situation in which 
S w {x w — X) 2 . This means that the maximum dis¬ 
tance (squared) between the electron inside the device 
active region and the surface is much smaller than the 
surface itself. In order to work out an approximate form 
for Eq. (31) in this regime it can be considered the fol¬ 
lowing change of variable x = ( x w — X). For simplicity, 
we assume that the electron is located on the left of the 
surface (i.e. X < x w —> x > 0) then: 


1 Jdp w p w V{p w f)x s ) _ J(x s ,t m ) 


V{x s ) 


\ip(x s ,t m )\ 2 


— n(x s , t m ), 


which is the result in Eq. (14) of the letter that we 
wanted to demonstrate. 


4>(\) = — tan 


, 4 X\/ 1 + 


2x 2 


Then, calling £ 2 = Eq. (32) becomes 


(32) 


B - DISPLACEMENT CURRENT AND 
MOMENTUM MEASUREMENTS ON A LARGE 
SURFACE 

In the text, we use the quantum version of the Ramo- 
Shockley-Pellegrini theorems in Eq. (8) to demonstrate 
that the measurement of the displacement current in a 
large surface is equal to the momentum measurement. 
Hereafter, we provide an alternative demonstration of Eq. 


$(£) = — tan 1 ( — | , 

\2y/e(l+?)J 


(33) 


such that the condition S w % 2 becomes equivalent 
to £ <C 1. So Eq. (33) becomes simply: 


*<««■«'=£ t-- 1 y?)- 


(34) 


















Remembering that tan 1 (a^)+tan x (^) = f for £ > 
0 then one has: 


$(0 


q r 7r 
7re . 2 


tan 1 (2£) 


(35) 


In Eq. (35) the term tan 1 (2^) can be expanded ob¬ 
taining: 


where v Xi \/i is the x-component of the Bohmian veloc¬ 
ity of the i-particle. One can reasonably assume that, 
for a large number of particles in the cable, we get 
(V4>(X fc (t)) • Vfc) « 0. So, finally we arrive to the 
result in Eq. (9), 


(It)s w oc (p Xl ). (40) 


<*>(0 


q 7T 

7T 6 2 


2£ + 


(2g) 3 

3 


(36) 


This last expression, Eq. (36), can be truncated at first 
order of £ for our large surface. Thus recalling the original 
variables one arrives at: 


$(X) 


q 7T 
7T6 2 



A') 


oc X. 


(37) 


Eq. (37) is an important results, it demonstrates that 
the flux of the electric held generated by a particle in a 
very large surface is proportional to the position of the 
electron. 

Now it can be discussed the general problem considered 
here, i.e. derive a microscopic analysis of the measure¬ 
ment of the total electrical current in a large metallic 
surface. In order to do that one has to “enlarge” the 
system considering also all the electrons composing the 
metallic surface, as done described in the main text. 
Without assuming nothing about the dynamics of the 
electrons in the metal, one can say that they contribute 
to the flux of the total electric held as described by Eq. 
(30) by superposition principle. One obtains, suppress¬ 
ing the dependence on x w and making reference to the 
position of the electron in the device as Ad, the following 
expression: 


N 

*(x 1} X 2 ,..., Xjv) =Ql!+a^ $(X fc ), (38) 

k—2 

where the actual Bohmian positions of the particles X,; 
have been used and a is a suitable constant. In Eq. (38) 
one can clearly see that the total electric flux is due to a 
contribution from the electron in the system oc Xi and 
another due to all the other electrons in the metal. 

So far, it has been considered that the electron in the 
active region of the device is not crossing the surface and 
therefore one gets that the total electric current is due 
only to the displacement current contribution. So the 
total current becomes: 


d ( A 

J Ts w oc ~rr = J a^W+a^$(Xi(f)) 

1 \ k—2 

N 

OC v Xl +^V$(X fc (t)) • Vfc, (39) 


Eq. (40) shows that the mean value of the total electri¬ 
cal current in a large metallic surface is proportional to 
the mean value of the momentum (x-component, i.e. the 
component perpendicular to the surface) of the quantum 
particle in the device. Let us emphasize that we arrive 
to the same result, Eq. (9) in the letter, from completely 
different arguments, without using the Ramo-Shockley- 
Pcllegrini expressions, nor the quantum expression of the 
particle current density. 


C -FREQUENCY DEPENDENCE OF THE WEAK 
MEASUREMENT OF THE MOMENTUM 

Here it is addressed the question on how the frequency 
on which the measurement is performed influenced the 
weakness of the measurement. In general a weak mea¬ 
surement requires that: the measurement does not per¬ 
turb (too much) the wave function, that the information 
extracted from a single experiment is not reliable and 
that the mean value corresponds exactly to the mean 
value of a strong measurement of the same quantity. 
In Fig. 2a in the letter it has been shown that if one 
wants precise information about the measured system, 
the perturbation of the wave function is increased, and 
vice versa. It has been also proven that the mean value 
obtain from the simulations with ammeter is equal to 
the mean value without ammeter (of the system alone), 
as clearly shown in Fig. 2. 

In the letter, we show numerically and explicitly that 
the weak measurement of the total current can be written 
in the language of Gaussian measurement Kraus operator 
as: 


^ r _ c p—p-w ) 2 

Iw(Pw) = c w dpe \p)(p\, (41) 

where p is the momentum (x-component) of the 
particle in the device, which is exactly the form used in 
the main text. 

Now, it will be discussed how the width <j w of the 
Gaussian measurement Kraus operator changes with the 
frequency of the current measurement. In Fig. 4 it is re¬ 
ported how a w varies with the frequency of the measure¬ 
ment. It can be seen that lowering the frequency yields 
to a more precise (mean) information about the system 
measured. As said in Ref. [27] our ammeter have to work 
at a frequency / < 1/r where r is the dwell time of the 
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FIG. 4. Blue dashed line probability distribution of the mea¬ 
sured total current at a frequency of / = 500 ■ 10 12 Hz. 
Red solid line probability distribution of the measured total 
current at / = 50 • 10 12 Hz. 

electron in the specific device, otherwise we will not get 
information about the interference pattern of the wave 
function in the numerical experiment considered in the 
main text. 


D - TOTAL CURRENT AND POSITION 
MEASUREMENTS ON A SMALL SURFACE 

In the letter, we have demonstrated through the use 
of the quantum version of the Ramo-Shockley-Pellegrini 
theorem that, in the case of a surface S s , where the elec¬ 
trical flux is collected is very small, compared to the 
distance of the electron from the surface, the measure¬ 
ment of the displacement current can be interpreted as 
a position measurement. Here we provide an alternative 
demonstration of Eq. (12) in the letter through the use 
of the Conditional wave function algorithm described in 
the letter. 

We consider the case in which the surface S s , where 
the electrical flux is collected is very small compared to 
the distance of the electron from the surface \ = x s — X. 
Specifically one has to consider Eq. (32) in the limit in 
which S s -C x assuming as before that the electron is 
on the left of the surface, i.e. x > 0- Making the fol¬ 
lowing change of variable £ = the previous condition 
becomes equivalent to study the condition £ <C 1 for the 
function 

(42) 

it can be easily seen that the first order expansion of 
the last expression is: 

= f “ *<*- < 43) 


The physical interpretation of this term is quite natu¬ 
ral: the contribution of a particle to the total flux mea¬ 
sured in a small surface is only relevant when the particle 
is “near” the surface. This means that the contribution 
is relevant when the electron cross the surface, i.e: 

X ss x s + \fSs with S s -A 0. (44) 

We note that the function (31) has a discontinuity in 
the point X = x s of the first kind. This is just an arti¬ 
ficial result due to the consideration of the displacement 
current alone. If we add the particle current in the dis¬ 
cussion of the total flux of the electric field, then the total 
current has no discontinuity. The particle current is just 
a delta function centered at x s , i.e. q5(X — x s ). 

Thus, since the flux of the displacement current can be 
interpreted (roughly) as a type of Heaviside step func¬ 
tion, its time derivative will be also proportional to a 
delta function. Therefore, the total (quantum ensemble) 
current measured on the small surface S s will be : 

{It) s a oc (6{xi - x s )) = 

dxiip*(xi,t)6(xi - x s )ijj(xi,t), (45) 

where ijj{x\ 1 t) is the wave function of the electron in 
the device. Eq. (45) explains the obvious relation that 
measuring the total current in a small surface provides 
information whether or not the particle is passed through 
the surface. In principle a treatment including the rest 
of the electrons of the system can be provided but in 
this case the result obtained does not change; in fact the 
electron in the device contributes to the total current only 
when pass through the small surface and then when the 
interaction with all the others electrons in the metal is 
strong (the distance between the electron in the quantum 
system and the metal electrons is very small, implying a 
strong perturbation ). 

Thus it is possible to write the measurement of the 
total current in a small surface S s in the language of the 
Gaussian measurement Kraus operator as: 

I s (x s ) = C s J dxe 2 \x)(x\, (46) 

which is exactly the form used in the Eq. (12) in the 
letter. 

E - CLASSICAL TREATMENT OF THE 
ELECTRIC FIELD 

The treatment of the electromagnetic interaction be¬ 
tween electrons is treated classically in the sense that no 
quantization of photons is considered. Although it is a 
typical assumption in most of the high frequency simula¬ 
tions of electron devices, hereafter we provide a detailed 
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justification. Following Ref. [31] we have that an electric 
field can be treated classically if the following condition 
is satisfied: 


\E\ 2 » 


h 

c 3 At 4 eo 


(47) 


This condition arise from counting the number of pho¬ 
tons in a time interval At. In our simulations we have 
that the time interval is the inverse of the frequency re¬ 
ported in Fig. 2, so At = 2 • 10” 14 s. For this value of 
the interval of time At the right hand side of Eq. (47) is: 


h 

c 3 At 4 eo 


2,7- 10 6 


N 2 

C 2 ' 


(48) 


We do not have, from the simulations, the direct infor¬ 
mation about the modulus of the electric field, but it can 
be estimated from the flux of the electric field through 
the surface used for detecting the displacement current. 
In particular from the simulations reported in the main 


text we have: 


\E\ = (49) 

D w 

where \E\ is the modulus of the mean electric field 
over the entire surface. In the simulations reported in 
the main text we have used a surface of S w = 10” 13 m 2 . 
The flux of the electric field, through the same surface, 
is $ « 7 • 10” 10 V ■ to. Thus for these values of the 
parameter we have 


\E\ 2 « 5 ■ 10 7 (50) 

satisfying the relation reported in Eq. (47) and thus 
justifying the use of a classical electrical field in our sim¬ 
ulations. We remind that, due to the dimensions of our 
device, we can use the (non-local) Coulomb potential as 
a reliable approximation in our simulations. This is a 
standard approximation for transport in mesoscopic sys¬ 
tems. 




